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ABSTRACT: A simple - yet plausible - model for B-type vortex breakdown flows is postulated; one that is based on the 
immersion of a pair of slender coaxial vortex rings in a swirUng flow of an ideal fluid rotating around the axis of symmetry 
of the rings. It is shown that this model exhibits in the advection of passive fluid particles (kinematics) just about all of 
the characteristics that have been observed in what is now a substantial body of published research on the phenomenon of 
vortex breakdown. Moreover, it is demonstrated how the very nature of the fluid dynamics in axisymmetric breakdown 
flows can be predicted and controlled by the choice of the initial ring configurations and their vortex strengths. The 
dynamic intricacies produced by the two ring + swirl model are illustrated with several numerical experiments. 
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1 introduction 



Vortex breakdown is officially five decades old this year: Although there is evidence to suggest that this phenomenon was 
observed several centuries ago in various atmospheric contexts, the first officially recorded description appears to be that in 
the experimental paper of Peckham & Atkinson P31 . This was followed hard upon by the more vortex breakdown focused 
investigations of EUe ll20l and Lambourne & Bryer ll37l . and what appears to be the first real attempt at analyzing the 
phenomenon in Squire 1561 . And even after half a century of research, interest in vortex breakdown has not diminished; 
in fact, it has probably increased significantly over the last fifteen or so years. 

The early work on vortex breakdown quickly captured the attention of a substantial group of talented fluid mechani- 
cians attracted by the importance, and challenge of unraveling the secrets of this intriguing and rather mysterious phe- 
nomenon. With this intense level of scrutiny, it did not take long to realize from the experimental evidence that there were 
apparently two distinct types of vortex breakdown configurations: an ostensibly axisymmetric bubble type (B-type) man- 
ifestation; and a definitely asymmetric spiral type (S-type) structure often observed as emanating from the trailing edge 
of a B-type breakdown form. Moreover, inspired by Squire's attempt to formulate a theoretical explanation of the phe- 
nomena in terms of variationally based critical state criteria, many similar and quite different theoretical descriptions were 
propounded. The three main types of theoretical characterizations employed critical state concepts, boundary layer sepa- 
ration analogies, and long wave hydrodynamic instability formulations, and a variety of flow properties and characteristic 
combinations of parameters, such as the Rossby number, were identified as playing key roles in the formation of vortex 
breakdown configurations. Outstanding examples of these contributions, in addition to Squire's pioneering analysis, can 
be found in the work of Benjamin Ludwieg Eo), Hall Leibovich lH, Escudier & Keller lEB, Trigub lEl, Spall 
et al. Il54l . Berger Q, Rusak et al. Il46l . Gelfgat et al. Il23l . and Krause 1341 [35l . These and other theoretical investigations 
have done a great deal to expand our understanding of vortex breakdown, especially the B-type variety. Nevertheless, the 
venerable enigma that is vortex breakdown has yet to yield to a complete and universally accepted exposition. 

Soon after the quest for a theoretical explanation of vortex breakdown began, a host of efforts directed at identify- 
ing vortex breakdown flows via computational techniques applied to both the Navier-Stokes and Euler equations were 
undertaken. Early results in this vein were somewhat inconclusive, but improvements in numerical methods and major 
advances in computers have produced results that are quite compelling. Notable examples of this important aspect of 
vortex breakdown research can be found in Grabowski & Berger l24l . Krause et al. l32]| . Hafez et al. [26 1, Spall & Gatski 
l55l . Breuer ITSl . Weimer fi601 . and Krause l33l . All of these numerical investigations have benefitted from and been 
complemented by a number of quite sophisticated experimental studies of vortex breakdown that began in earnest around 
the same time, with Sarpkaya l48l . Faler & Leibovich l22l . and Uchida et al. \59\ among the best examples of the earlier 
experimentally based investigations. 

Holmes l29l was instrumental in starting a relatively new trend in research in fluid mechanics in general and vortex 
breakdown flows in particular: During the last twenty years, many investigations of a theoretical, computational, or 
experimental nature have been conducted from the perspective of modern dynamical systems theory. Terminology such 
as structural stability, stable and unstable manifolds, homoclinic and heteroclinic orbits, and strange attractors have now 
become the lingua franca for a large segment of vortex breakdown investigations. Prime examples of this trend can be 
found in the work of Blackmore iH, Blackmore & Knio lT2l . Br0ns and his collaborators lfT6l [TTl [Tsl . Gelfgat et al. 
l23l . Serre &. Bontoux l49l . S0rensen &. Christiansen ISTll . and Sotiropoulos and his coinvestigators |l52]|53l. Viewing the 
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results of this dynamical systems approach, it is safe to say that it has already proved quite useful in apprehending some of 
the more complex features of vortex breakdown flows, and has great potential for catalyzing future leaps of understanding 
in this area. 

Dynamical systems theory certainly inspired the work presented here on our vortex breakdown model, which was first 
adumbrated in ITZi . But the idea of using a pair of circular rings was also influenced and guided by the experimental and 
computational literature and interaction with some of the leading practitioners in these areas. On the one hand, we knew 
from the dynamical systems perspective that, among other things, two rings are capable of generating the kinds of chaotic 
kinematics that appear to be quite common in vortex breakdown flows lH \lU\ (TT] [T3] [SO], but on the other hand, our 
choice was informed by the apparent ubiquity of two coaxial vortex rings in numerical and experimental investigations. 
Naturally, our decision to concentrate on an ideal fluid based model was motivated by a desire to keep the model simple 
and also have the consequent rich symplectic structure at our disposal; the advantages of which are, we believe, amply 
demonstrated in the sequel. 

The development and analysis of our two coaxial vortex ring model for B-type vortex breakdown unfolds in this 
paper as follows: In Section 2, we formulate the mathematical model for the motion of a pair of coaxial vortex rings 
immersed in a swirling ideal fluid flow along the axis of symmetry of the rings in the context of Hamiltonian dynamics, 
wherein we consider such features as the complete integrability of the resulting system. Next, in Section 3, we derive the 
Hamiltonian dynamical equations for the motion induced in passive fluid particles (kinematics) by the dynamics of the 
rings in the swirling flow. The associated kinematic equations are formulated in two ways: in terms of a passive third 
ring with zero vortex strength; and directly as a (time-dependent) Hamiltonian one-degree-of-freedom system. It is noted 
that the resulting kinematics is not completely integrable, so it is capable of exhibiting chaotic flow regimes. In addition, 
we introduce the approach to studying the dynamics via a Poincare map associated to successive intersections of the 
streamlines with a fixed meridian half-plane. Then the usual aspects of the dynamical system for the ring motion, including 
classification of stationary points, are analyzed in Section 4. This is followed in Section 5 with an analogous investigation 
of the kinematic equations, which incorporates a brief analysis of homochnic and heterocUnic orbits connecting the several 
fixed points. 

Since it was found that the kinematics for the case of fixed rings is too regular to subsume the kinds of exotic streamhne 
patterns that one expects in vortex breakdown on the basis of several careful experimental and numerical investigations, 
we formulate, in Section 6, the model equations in terms of small (quasiperiodic but not necessarily periodic) oscillations 
of the rings. Particular attention is paid to obtaining asymptotic expansions of the heterochnic orbit comprising the outer 
boundary of the bubble for the trace of the kinematics in a meridian half-plane. Then, in Section 7, we prove using 
Melnikov's method that small oscillations in the rings produce chaotic streamline configurations of the type observed - 
but not proved - in earlier studies of vortex breakdown phenomena. This is followed by Section 8 in which a variety of 
numerical experiments are run to illustrate the dependence of the kinematics on various model parameters such as the 
swirl strength and relative vortex ring strengths. The nature of transitions to chaos as the size of the ring oscillations is 
increased is given special emphasis. Finally, we conclude in Section 9 with a discussion of our model, its plausibility, 
and its efficacy in predicting and controlling vortex breakdown behavior We also indicate some natural further research 
directions suggested by the results obtained in our investigation. 
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2 The Two Ring Model 



First, we consider the dynamics of two coaxial vortex moving in an ideal (= inviscid and incompressible) fluid initially at 
rest in K'^. The axis of symmetry of the vortex rings is chosen to be the x-axis, and their respective positive strengths are 
denoted as Fi and F2. Owing to the axisymmetry of the motion of these rings, their positions are completely determined 
by their intersection points with any meridian half-plane containing the x-axis. In order to fix ideas, we choose this half- 
plane, which we denote by S), to be the xy-plane with y > (j/ > 0, z = 0), so that the motion of the rings is characterized 
by their respective points of intersection :— {xi,yi,0) and {x2,y2) '■= {x2,y2,0) with the designated half-plane. 



Defining s 



+ z^, which reduces to on f), it is well known |0] [TO] HT] [T3] [161 [iTl that the equations 



of motion of the rings - desingularized to eliminate infinite self-induced velocities - may be expressed as the 2-degree-of- 
freedom Hamiltonian system 

"''/^ cos 2cFda 



KkSk = 4:KjKkrjrk {xk - Xj) 



KkXk 



for j,k — 1,2, with j ^ k, Kk 



2rk 

rfe/2^, 

A12 : 



log 



8rk 



3/2 



(1) 



''/^ {rj - Tk cos 2a) da 

' 



(''I - ^'2)^ + {xi - X2f' + 4rir2 sin^ cr, 



(2) 



< (5 ^ 1 is a very small positive number representing the common radius of the two rings in the desingularization 
procedure, and 

7 := i (^1 + log 2 + ^ e-« log ^d^^ ^ 0.558, (3) 

which is an artifact of assuming a Gaussian vorticity distribution in the cores {i.e. in the tubes of radius 5) of the rings 
in the desingularization. It should be remarked here that different desingularization approaches produce slightly different 
parameters and dynamical equations, but the overall qualitative aspects of the motion are unaffected by such choices, and 
the differences in the quantitative properties are quite small. 
By introducing the Poisson bracket 
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which is quite often employed in vortex dynamics, (1) can be recast in the manifestly Hamiltonian form 

Sk = K'f:^dj;^Ho = {Ho,Sk}, Xk = -n'^^dsf^HQ = {Ho,Xk}, {I < k < 2) 
where the Hamiltonian function for the system is 



(4) 



(5) 



Ho 




log 



8rfc 



(I + 7) 



4KiK2rir2 



7r/2 



cos 2ada 

A 1/2 



log(^)-(l + 7) 



+ 2kiK2 {ri2+ + ri2-) [F (A12) - E (A12)] 



(6) 
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where 

ri2± := (7-1 ± r-s) + (zi - za) , \i2 ■= \ , (7) 

and F and G are, respectively, the complete elliptic integral of the first and second kind given as 

r^^ d<7 r/'^ , 

F (A) := / = and E (A) := / VI - A2 sin^ ada. (8) 



This use of elliptic integrals in the equations of motion is classical and is often used {cf. lfT0l[30l[36l ). We note here that 
(5)-(6) is (real) analytic on the following subset of M^, which serves as the phase space of the system: 

Xd ■= {isi,S2,Xi,X2) e : si,S2 > 0, (si,xi) ^ (s2,a;2)} • 

The autonomous Hamiltonian system (5)-(6) has the following independent integrals in involution (see e.g. ll2ll4l[T0l 
[12 mi HI): 

2 2 

ifo, =;^Kfcr2, (9) 

k=l k=l 

SO it is completely integrable in the sense of Liouville-Arnold (LA-integmble). We note here that a system of two or 
more coaxial rings - even if one ring has a vortex strength of zero - does not have any additional independent constants 
of motion in involution, as proved by Bagrets & Bagrets Consequently, the dynamics of three or more rings is apt to 
include chaotic regimes in most cases. 

2.1 Full model including swirl 

We now immerse the two ring model in an ambient swirling flow. To be more precise, we shall assume that the rings start 
and remain in a swirling flow about the x-axis. In order to insure that the dynamical equations for the fuH model including 
the coaxial vortex ring pair and the ambient flow still have a Hamiltonian formulation, we shall choose our models for the 
ambient swirling flow to be axisymmetric potential flows. A convenient way to represent such swirling flows is via the 
use of cylindrical coordinates having the x-axis (rather than the usual z-axis) as the axis of symmetry. These coordinates 
naturally are defined as 

a; = X, y ~ rcos9, z — rsin^, (10) 

where the angle 9 is measured counterclockwise in the yz-plane starting at zero along the positive y-axis. The general 
expression for the ambient velocity in terms of these coordinates is 

Va = (ia,ya,2a) = {Ua,Va,Wa) = (^Xa, r a COS 9 - ZaO , ra Sm 9 + yaO^ ■ (11) 

We assume for the ambient swirling flow that Xa is a quadratic function of s = of the form 

Ua (s) = Xa (s) = -a (l + fliS + 025^) , (12) 

where a > and ai and 02 are real constants to be chosen in the sequel. To preserve axisymmetry, we assume that 9 is a 
polynomial function of s = of degree two. Whence, we see that the ambient velocity assumes the form 

Va = {Xa,ya,Za) = (^Xa, Ta COS 9 ~ Za9 , ra Sin9 + ya9j = (^-Q (l + ai.S + 025^) , ~ Za9{s) , ya9{s)^ , (13) 
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where we shall only consider the following polynomial forms for the angular velocity 

e {s) = n (1 + bis + b2s^) , 



(14) 



where there is no loss of generality in assuming that > 0, and bi and 62 are real constants to be specified in the course 
of our investigation. 

Now we are in a position to formulate the full coaxial vortex ring pair + swirling flow dynamics. Once again, the 

axisymmetry enables us to specify the motion, modulo rotation about the x-axis, of the rings in terms of their intersection 

points with the meridian half-plane Sj. Observe that the "modulo a rotation about the x-axis" is not really a problem in 

describing the evolution of points on a ring. To see this, suppose we consider a point on a ring, say the first ring, that is 

initially at (a;i(0), si(0), 0i(O)), and we know xi{t) and si{t). Then we can find 9i{t) from (14) and a simple integration 

as ^ 

91 (t) := e (si (t)) = 01 (0) + Q f (1 + ciSi (r) + c^sj (r)) dr, (15) 

Jo 

and there is an analogous formula for 62{t) := 6{s2{t)). 

Adding the ambient swirl to ( 1 ) only changes the equations of motion by a the addition of a simple polynomial function 
along the x-axes; namely, the dynamical equations become 



{Xk -Xj) 

Jo 



^^"^ COS 2(Tda 



A 



3/2 
12 



KkXk = -Kka (1 + aiSk + a2sl) + 
having the following Hamiltonian analog of (5)-(6) 



log ( ) - 7 



+ 2k 



jKkVj 

Jo 



^^"^ {vj — Tfe COS 2cr) da 
^12 



(16) 



Sk = Kf, ^dx^H = {H,Sk} , Xk = -Kfe ^ds^H = {H,Xk} , (1 < < 2) 



with Hamiltonian function 

2 



ai 2 ,02 3 



H := a^Kk [sk . 2 ■ 3 

fe=i 

' 2 

Y^'^l'Tk 



,2 , "2 ^3 



fc=l 



log(^)-(l+7) 



+ 4Ki«;2rir2 



''/^ cos 2tT(icr 



Kk [Sk + ySfe + ys; 



(17) 



k=l 



, fe=l 



+ 2K1K2 (ri2+ + ri2-) [F (A12) - E (A12)] } . 



(18) 



Observe that (17)-(18) also is LA-integrable inasmuch as it has the following pair of independent invariants in involution: 



H, G := ^ KkSk = ^ Kkvl- (19) 

k=l fe=l 

Our primary concern in the sequel is the affect of the ring dynamics - as described by (17)-(18) - on the motion of 
passive fluid particles, which is often described as the fluid kinematics. 
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3 Dynamics of Model Kinematics 



The dynamics of passive fluid particles induced by the motion of the coaxial vortex ring pair if the swirling flow charac- 
terized by (17)-(18) can be described in two primary ways: Either indirectly and autonomously in terms of a restricted 
three ring problem in which we treat a passive fluid particle as belonging to a third ring of zero strength, or directly as 
a time-dependent system of ordinary differential equations in which we first solve for the dynamics of (16) in order to 
describe its affect on a passive fluid particle. Each of these approaches has its advantages, and we shall find it convenient 
to develop both of them in what follows. 

3.1 Passive third ring approach to the kinematics 

In this approach, we actually consider a three ring problem with two of the rings as above with their given nonzero 
strengths, Fi and and a third (advected) ring of strength F = 0. Once again, the axisymmetry implies that the 
motion of the rings is completely determined by that of their respective intersection points, (xi.yi), (x2, ^2) and (x, y) 
with the half-plane Sj. The equations of motion of these three points are analogous to (16), and can be written as(c/ 

iiiioiiniiniiisiEi) 



, cos2ada 

KkSk = 4:KjKkrjrk(Xk - Xj) J72~' 

Jo 



Kfcifc = -Kka (1 + aiSk + a2sl) + ( 



s = r 



X Xk J 



fe=l 



cos 2ada 



,3/2 



X ~ —a (1 + ais + a2S^) + 2 



Kkrk 



log ( ^ ) - 7 



fe=i 



7r/2 



(rj — rk cos 2cr) dcr 
2KjKkrj I , 

^12 



{Tk ^ f cos 2(t) da 



(20) 



where 



Afc := (r - rk) + (x - Xk) + Arvk sin^ cr. (1 < fc < 2) 



(21) 



This system can be recast in what amounts to a kind of piecewise Hamiltonian form. In particular, we may rewrite 
(20) as 



Sk^ K]:^d^^,H = {H,sk} , Xk = -K^^ds,H,^{H,Xk}, (1 < A: < 2) 
s = d^H^ = {H^,s}^ , X = -dsH^ = {H^,x}^ , (22) 



where {•, •} is as in (4), H is as in (18), and 
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H^:=H + a 



4^E 



Kfcrfc 



k=l 



ir/2 



COS 2(jda 

A 1/2 



fc=i 



4KiK2rir2 



4r > ^ Kfcr-fc 



1/2 



fc=l 



fc=l 
""/^ COS 2ada 



log(^^)-(l + 7) 



1/2 



ai o 02 ^ 



.2 + 21.3 



, fe=i 



log 



S 



(1 + 7) 



2kiK2 (ri2+ + ri2_) [i^ (A12) - (A12)] + 2 ^ Kfe + rfc_) [F (Afc) - £; (A;, 



fc=i 



where 



:= (^fc ± rf + (z,. - zf , A^ := '^'+ , , (1 < fc < 2) 



(23) 



rk+ + rfe 

F and E are as in (8), and the new Poisson bracket is actually the standard one defined as 

^^'3^*- \dxds dsdxj' 
The system (22)-(23) is analytic on the six-dimensional symplectic space 

Xdk := {{si,S2,xi,X2,s,x) £ M*^ : si,S2,s > 0, (si, S2, xi, 2:2) G X_d, {s,x) ^ {sk,Xk), 1 < fc < 2} , 

and clearly shares the following two independent constants of motion with (16): 



(24) 



(25) 



(26) 



fc=i 



fc=i 



but no additional integrals in involution as demonstrated in [|4J. From this we infer that the motion of the advected ring 
can become chaotic, which is consistent with the apparent chaotic motion found in vortex bubbles llT6l[T5l[T7ll23l[39ll49l 



3.2 Direct formulation of kinematics 

In the direct method, we obtain a time-dependent system of differential equations for the motion of a general passive fluid 
particle by assuming that the dynamics of the pair of coaxial rings has already been determined; i.e., we have the solution 

{si{t),xi{t)), {s2{t), X2{t)) of (17)-(18), which means that we also know 6i{t) := 9{si{t)) and 02{t) := 0{s2{t)) owing 
to (15). 

We could start with the Biot-Savart law as was done in Blackmore & Knio 1 10] and derive the system of differential 
equations governing the motion of points x = (x, z) G M'^. However, we simplify things by taking advantage of the 
fact that the rotational motion of the particle about the x-axis is completely determined via (14)-(16) by the motion of the 
intersection point with of the advected circular ring containing the particle. Not only does this reduce the dimension of 
the (continuous) dynamical system from three to two, but it also has the advantage of having a Hamiltonian formulation. 
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albeit a time-dependent one. Taking this approach, it is easy to see that the desired system is akeady at hand - namely 
from the last two equations of (20) with the positions of the two rings specified as functions of time. In particular, the 
equations of motion of the intersection points of the advected ring with ^ are 



2 

s = 4r KkTk {x - Xk (i)) / 



"372 



, , kTkit) f 



x = -a{l + ais + a2S^)+2Y,Kkrk{t) 1^ ^ ^ ^ ^ ^ — , (27) 



where, 

Ak{t):={r-rk{t)f + {x-Xk{t)f + 4.rrk{t)sm^a. {l<k<2) (28) 
The system (27) is (time-dependent) Hamiltonian since it can be expressed in the form 

s = dxH {s, x,t) , X = —dsH (s, x, t) , (29) 

where the Hamiltonian function is 

rV2 j,Qy2crrfa 



n {s, X, t) := a (s + ^ ^^sA _ 4^ V KkVk (t) I 
^ ^ ^ fe=l -'o 



it) 



= a(^s+^s^ + ^ 2 ^ /^fe {rk+ (i) + Vk- {t)) [F (Afc(t)) - E (Afe(O)] , (30) 

and we have used the notation in (24) interpreted as (known) functions of time. We note that this system is analytic on the 
phase space 

Xk {{s,x,t) e : s > 0, is,x) ^ {sk{t), Xk{t)) ,l<k<2}. 

Observe that it is clear from either (20) or (27) that if the passive particle starts on the .T-axis (s = 0), it remains there for 
all time. To be more precise, the hypersurface s = is an invariant hypersurface of the six-dimensional phase space for 
(20), and it is an invariant line in the trajectory space associated to (27). From a fluid dynamics perspective, the points on 
s = for which a; = are of particular importance since they correspond to stagnation points on the axis of symmetry, 
which one would expect to define the leading and traiUng edge (points) of a B-type vortex breakdown structure. We shall 
have much more to say about this in the course of our analysis, which follows. Another observation that will prove useful 
for the analysis, is that the half-plane plays a very natural role as a transversal for a fairly obvious Poincare map - a map 
that conveys a great deal of useful dynamical information. As a preview, we give a brief description of this map. Let 
{s{t),x{t)) be the uruque solution of (27) starting at {so,xo) € Sj, with sq > 0. Then the trajectory in winds around 
the X-axis until at a first ti > 0, it intersects Sj again. This establishes a smooth mapping of Sj into itself via the formula 
n {{sq, xq)) := {s{ti),x{ti)). This handy mapping shall be investigated further in the sequel. 



4 Analysis of the Ring Dynamics 

As for the ring dynamics, governed by (17)-(18), it was pointed out in the preceding section that it is LA-integrable, a 
property often referred to as integrable by quadratures. Therefore, at least in a theoretical sense, we can integrate the 
dynamical equations (17)-(18). It appears, however, that this integration cannot generally be carried out in a practical way 
in order to find a tractable closed form solution for the equations of motion, except perhaps in some very special cases. 
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4.1 Stationary points 

We begin with an investigation of the fixed points of the system, which we consider in the form of (16). Solutions of 
si = S2 = are easy to characterize. This follows from the readily verifiable fact that the integral 



Jo 



^12 



(31) 



is positive for all choices of ri, r2 > 0, xi and X2 ( except in a limiting case where one ring collapses to a point with its 
core radius decreasing in a special way). Thus a necessary and sufficient condition for = S2 = is that xi = x^. This 
connmon value shall be denoted as ^ := a;i = X2. 

Naturally, it remains to find solutions of xi= ±2 = 0, and to this end we introduce some convenient simphfications. 
There is no loss of generality in assuming 



1 = Ki < K := K2 and < ri < r2 



(32) 



since this can be attained by a rescahng of the time variable t and reordering the variables, if necessary. In order to 
simphfy matters, we shall assume in (12) that ai= a2 = 0, so the equations we need to solve take the form 

'""'"^ {r2 -ricos2a)da 



1 



= ±1 = -a+ — 
2ri 



= X2 = -a + - — 
2r2 



log(^ 



■7 



2Kr2 



^3/2 
^12 



log ( ^ ) - 7 



2ri 



.77/2 

Jo 



{i^i — T2 COS 2(t) da 



,3/2 
^12 



where 



Ai2 := (ri - r2) + 4rir2 sin^ u. 



(33) 



(34) 



It is convenient to display the full set of equations defining the fixed points, which owing to (16) and our assumptions can 
be written as 

t/2 



^ - / ^ r' cos2ada „ 

si = ^i{ri,r2,xi,X2;a,K) := 'lKrir2{xi - X2) / ?75— = 0' 

Jo A?/' 



^ - / N cos2ada „ 

S2 = ^2{ri,r2,Xi,X2;a,K) := 4:rir2{x2 - xi) / 372"^^' 

Ai 



'12 



1 



xi = *i {ri,r2,xi,X2;a,K) ■= -a + — [log(xri) -7] +2Kr2 



L 



''''^ (r-2-ri cos2a)d<T 
^12 



.3/2 



= 0, 



K f 

X2 = *2 {ri,r2,xi,X2;a,K) :=-a+ — [log(xr2) - 7] + 2ri / 

2r2 Jo 



""^^ in -r2 cos 2a) da 

^12 



= 0, 



(35) 



where % ;:|> 1 is a multiple of the curvature of circular cross-sections of the boundary of the (virtual) core of each coaxial 
vortex ring, and is defined as 



(36) 



We note that, strictly speaking, the functions in (35) also depend on the parameter x; however, given the "virtual nature" 
of the definition of x> we choose to fix its value at 1,000 and suppress it as a coordinate of parameter space in the 
sequel. As the first two equations of this system require that Xi = X2 ■= ^, this solution reduces to the following pair of 



10 



^-independent radii equations: 

1 r / X 1 f^^^ ('2 - ''1 (-ori ^(7) da 
*i {ri,r2;a,K) :=-«+ — [log(xri) - 7J + 2Kr2 / ^^y^ = 0, 

* / X n / N 1 „ /■'^^^ (ri — r2 COS 2o-) da „ 
*2(ri,r2;a,K) :=-a+ — [log(xr2)-7] + 2riy^ l^f"^ " 

4.1.1 Solution of radii equations 

It can be shown via a straightforward argument, in which the function 

V(r) :=(2r)-Mlog(xr)-7] (38) 

plays a key role, that the radii equations (37) have precisely four solutions (ri,r2) when a lies between the value one 
and a not too large upper bound. In these four solutions, there is one with both coordinates very small, one with both 
coordinates significantly larger than in the first solution, and two solutions with one coordinate small and the other quite 
a bit larger. We leave the proof of this to the reader, and only give a sketch of the reasoning involved. Toward this end, it 
is useful to note some basic properties of the graph of the function ^ restricted naturally to r > 0. It is easy to see that 
ip has a unique zero at ro := e'^x~^ its derivative ip' also has a unique zero at r = p := e^'^~^^^x~^ which the 
function attains its maximum of tp (p) = 2~^e~^^~^^^x = P~^- Moreover, ip' (r) is positive for < r < p and negative 
for p < r < 00, and il>{r) — > as r ^ 00. 

For our sketch of the solution of the radii equations, we shall find it instructive to ignore the integral terms in (37) 
thereby obtaining the reduced radii equations 



{ri,r2;a,K) :=-a+ — [log(xri) -7] = 0, 
*2 (n, ^2; a, k) := -a + ^ [log (xr2) - 7] = 0. (39) 

It turns out that the solutions of the full equations (37) are completely analogous to those of (39), and can be obtained 
rather directly using identical methods. Accordingly it follows from the form of (39) that we ought to restrict a in both 
(39) and (37) so that < a < p~^, which is not much of a restriction at all since p~^ is so large. It is convenient to recast 
(39) in the form 

- {2ari + 7) + log (xn) = 0, 
- (2ar2 + Kj)+K log (xr2) = 0. (40) 

It is easy to show that for 1 < a < p~^, each of (40) has a small (positive) solution and a larger solution to the left and 
right, respectively, of p. Let us denote these solutions, which are functions of (a, k), as fi^ , fi^ , and ^2i,> and note that 
is a simple matter to verify that they satisfy the following properties: 

ro < r^2„ < ri^ < p < fi, < f2(,, 

-r2<.(Q!,K),r26(a,K) -^^(q;,^) and r2fc(Q!, k) - ri„(a, /t) 
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are decreasing functions of a for every k > 1 and increasing functions of k for each a > 1, and finally 

f2a, («) k) i To as K t 00 

for each a > 1. We also note that it is easy to show that can be obtained as the limit of the following iterative scheme 
in which the first equation of (40) has been reformulated as a fixed point problem 

r("+i) = lexp(2ari"^+7), (41) 

where the initial iterate r^^ can be chosen as any value between zero and p. Moreover, ri_^ can be obtained as the limit 
of the following iterative scheme in which the first equation of (40) is reformulated as the (inverse) fixed point problem 




where the initial iterate r\ ' in this case can be chosen to be any reasonably large number greater than p\ for example, 
r^^ = 10 is always a suitable choice. Naturally, there are completely analogous schemes for finding f2„ and using 
the second of the equations (40). 

These observations about the reduced form (39) of (37) provide the basic tools for solving (37). More precisely, we 
find that although the additional integrals in (37) break some of the symmetry inherent in (39), there are all together four 
solutions, which we denote as (ri(i), 7=2(1)) = (?^i(i) (",«:), 7=2(1) (fi(ii) , 7=2(11)) = (^'^ii) (a, «) , ^2(11) (",/«)), 

(^i(iii),?^2(iii)) = (^i(iii) (a, k) , ^2(111) (a, k)), and (ri(iv) , r2(iv)) = (7^i(iv) (a, k) , f2(iv) (a, «)), and for which the 
following properties - somewhat mirroring those of the solutions of (39) - hold: 



< h{i) < h(i) 4fi(i); < f2(ii) -C fi(„) ~ (2 X 10^) f2(ii); 
< ri(iii) «; f2(iii) ~ (5 X 10^) fi(iii); < (lO"^) ri(i) ~ fi(iv) < hiw) - 
and f2(i) - ri(i), fi(ii) - hiii), ^2(111) - 'ri(iii), and hiw) - ^i(iv) 



(43) 



are all decreasing functions of a for every k > 1 and increasing functions of k for each a > 1. We also note that all four 
solutions can be obtained iteratively using the procedures described in (41) and (42), although considerably more care 
must be exercised with regard to selecting the initial point in the process. In particular, the solution (fi(iii) , 7=2(111)) may 
be obtained via Picard iteration in the form 




— < K log [xrf 



. (ri) (ri) 



^7r/2 (4"^ -rJ"^cos2cr) da 

7r/2 (^r^"' - r'^^'^ cos2a^ da' 
k%'{rt\rt') 



(44) 



with {rf\/^^^ = (0.001, 1.0), and then 

/ (n) (n)\ . N 

[r\ ',r2 j (ri(iii),r2(iii)j as n ^ 00. 
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4.2 Stability of stationary points 



We have now identified a whole 1 -parameter family of stationary points for (16) that can be represented as in terms of the 
parameter ^ as 



£ = ;= |g = {si,S2,Xi,X2) ■■ Si = si(i) ri(i),si = Si(ii) 



1(111). 



si = si(iv) := ^i(iv)' = S2(i) := ^2(1)' ^2 = S2(ii) := f^u), S2 = S2(iii) := ?^2(iii). (45) 
S2 = S2(iv) := rl(iyj,xi =X2=^, ^ G m| . 

To perform a standard linear stabiUty analysis for a typical e £, we compute the derivative (matrix) for X12 := 
($1, $2, *i, *2); namely, 

a,, $2(9) a,, $2(9) a,i*2(g) a,. $2(9) 

ds,^i{q) d,,^i{q) 
Va«,*2(9) a«,*2(g) d,,^2{q) d^,^2{q) ) 



^'12 (9) 



(46) 



In aid of this, we compute that 



(q) = 0, ds.^i (q) = 0, (g) = 4«fin 



7r/2 



ds,^2 (q) = 0, a«,$2 (g) = 0, a,, $2 (g) = -4fif2 / 

JO 

(g) = [log (xn) - (1 + 7)] - 



cos 2(7 da" 

A%^{fi,f2) 

""l"^ cos2ada 



9x2^1 {q) = -4«;rif2 



''/^ cos 2CTdcr 



9x2^2 (g) = 4fif2 / 

Jo 



Aijf(f 1,7=2)' 

''/^ cos 2CTdcr 



A?f(fi,f2) 



2Kf; 



^72 ^Ai2 (ri, f2) cos2cr + 3 (f2 — ri cos2cr) [(fi — f2) + 2r2 sin^ cr]^ dcr 



^ /•7r/2 [(1 + f2 - n cos2it) A12 (ri,r2) + 3 (r2 - fi cos2cr) [(r2 - h) + 2fi singer] j rfcr 



(fl,f2) 



5x,*i(g) = ax,*i (g) = 0, 

2 |-w/2 ^(1 + fi - f2 cos2fT) Ai2 (ri,r2) + 3 (ri - f2 cos2cr) [(fi - r2) + 2f2 sin^ a]^ da 



ds,-^2{q) = ^ , 

ri Jo 

5x,M/2 (g)=5x,vi/2(g) = 0, 
2f2 \ 2r^ 



a.,*2 {q) = -. 



2 [l0g(xf2)-(l+7)] 



2ri 



7r/2 ^Ai2 (ri,f2) cos2cr + 3 (fi — f2 cos2(t) [(f2 — ri) + 2fi sin^ cr]^ tier 

(ri,r2) 



(47) 



To simplify our analysis of the spectrum of X'12 (g) at each of the fixed points, we introduce the notation 

Qi = QiiO ■= (si(i).S2(i),^,0 ' In = QiiiO ■= (si(ii),S2(ii),^,^) , 
9iii = qui (0 := (si(iii),S2(iii),^,0 ' ^^'^ 9iv = 9iv (0 := (si(iv)> S2(iv), C> • 
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It is a routine - albeit tedious - matter to show that all four of the derivatives X'12 {<li) , ^'12 {Qiv) have a pair of 
zero eigenvalues, which is consistent with the fact that none of the stationary points is isolated and all of the derivatives 
are symplectic matrices, and a purely imaginary conjugate pair: Asji) = ivi, A4(i) = A3(i) = —ii/i; X^u^ = iuu, 
^4(11) = A3(ii) = -ivii; A3(iii) = iuiu, A4(iii) = A3(iii) = -iuiu; and A3(iv) = ii^v, X^iy) = X^iw) = -iiAv, 
respectively, with vi,vii, vm, Piy > 0. 



5 Analysis of the Ring Kinematics 

With the ring dynamics disposed of in the preceding section, it is now quite simple to resolve the ring kinematics. We 
begin our analysis from the passive third ring approach as described in Subsection 3.1. Our calculations are rendered 
simpler owing to the fact that the first four equations of (20) do not depend on the last two equations, so the analysis of 
stationary points for the ring dynamics can be used directiy for the kinematics. 

5.1 Stationary points in the passive third ring approach 

First we rewrite (20) incorporating the assumptions of Section 4 and introducing some additional notation: 



si = ^i{ri,r2,xi,X2;a,K) := 4Krir2[xi - X2) / 372—, 

Jo 

''/^ cos 2(Tda 



S2 = ^2 {ri,r2,xi,X2;a,K.) := 4rir2 {x2 - Xi 



.3/2 
^12 



''/^ (r2-ri cos2cr)d£7 

372 ' 



3/2 
^12 



1 f 

xi = *i {ri,r2,xi,X2;a,K) :=-a+ — [log(xri) - 7] + 2Kr2 / 

^ri Jo 

T / X «^ n / X T « (''i ~ '^2 COS 2a) da 

a;2 = (ri,r2,a;i,a;2;a,K) := -Q! + — [log(xr2) -7J + 2ri / 

^r2 Jo 

^ / , ^ , , r^"^ cos 2ada 

s = ^{ri,r2,r,xi,X2,x;a,K) \= Ar y. i^kTk {x - Xk) / 

fc=i 

2 ^7r/2 

X = "1! {ri,r2,r,xi,X2,x;a,K) := -a + 2^K,krk / 

fc=l -^0 



,3/2 ' 



(rfe — r cos 2a) da 



(48) 



where ki — 1 < k := K2- To find the set of stationary points, we set all of these equations equal to zero. But we have 
already taken care of the first four equations; namely, we found the line £ defined by (45) comprised of a 1-parameter 
infinity of (stationary) points for which $1 = $2 = ^'i = *2 = 0. This leaves us to solve 



$ := 4r {x - ^) 
* = -a + 2 



ri 



'''^ coa2ada 



A 3/2 
A{ 



Kr2 



cos2ada 



a: 



3/2 



0, 



.77/2 

Jo 



(ri-rcos2a)da ^ r'^ 

'"0 Ap ^'^'Jo 



{r2 — r cos 2a) da 



= 0, 



(49) 



where 



(r 



fk)^ + ix-i)'^ +4rfksm^a, {l<k<2) 



(50) 
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for all choices of pairs of fi andr2 included in £; namely, ^^1(1),- • • , '^i(iv) andf2(i),. . ., r2(iv). respectively. It is therefore 
helpful to identify the following cases for the coaxial ring fixed point types: 

Type I. n = ri = fi{i), r2 = h = ^2(1) 

Type II. n = fi = ri(ii), r2 = r2 = ^2(11) 

Type III. n = fi = fi(i„), r2 = f2 = f 2(111) 

Type IV. n = fi = ri(iv), r2 = r2 = hiw) 

As the bracketed term in the first equation of (49) is positive, the only choices for zeros are r = and a; = ^. If r = 0, 
and we consider Type I, the second equation of (49) reduces to the form 



^ = -a + TT 



'i(i) 



Kr: 



2(1) 



fim + - 



3/2 



^i(i) + - 0' 



3/2 



(51) 



which immediately yields a pair of solutions that are symmetric with respect to any chosen value of a; = ^; namely, 



„(±) 



(52) 



where ?7i > 0. We note that ajj"*"^ and x\~^ correspond for Type I, respectively, to the stagnation points on the leading and 
trailing edges of the vortex breakdown bubble (on the axis of synometry). Analogously, we find points x^^^ := ^ ± rju, 
x[^ := ^± Tfiii and x^y^ := ^ ± rfiv for Types II, III and IV, respectively. 



If we choose a; = ^ for Type I (for a > 1) , the second equation of (49) becomes 



-a + 2S (r; f i(i) , f2(i) , = 0, 



(53) 



where 



S (r-;ri(i),r2(i),«;) := 



tt/2 



fl(I) { 



ri(i) 



r cos2ct) 



«7'2(i) [r2(i) 



r cos2cr) 



{r - h(i))^ + 4rfi(i) sin^ a 



3/2 



{r - ^2(1))^ + 4rf2(i) sin^ a 



3/2 



da. (54) 



Then a routine - but rather laborious - analysis of the properties of S shows that it is a smooth, nondecreasing function 
of r for ri(i) < r < r2(i) such that linirif^^j^ S = —00, limrff2(i) S = 00, and the derivative with respect to r increases 
rapidly from nearly zero to +00 very close to f2{i) ■ It also is easy to verify that H is nonvanishing when r ^ , '^2(1)]- 
Therefore, in light of the definitions of fi(i) and f2(i) as a solution pair of (39) with fi(i) < r2(i) — 47^i(i), we infer that 
(54) has a unique solution fj, which is increasingly near to fi(i) as a increases; in particular, 

< 'ri(i) < n < f2(i) and n J. fi(i) as a t 00. 

Observe that one should also include the singular stationary points (,s, x) = (si(i), C) and (s, x) = (s2(i), which are 
not actually included in the phase space - but are in fact stationary with respect to the kinematics - and can readily be 
shown to behave like (singular) centers surrounded locally by periodic orbits. It is straightforward to verify the analogous 
behavior for solutions of (54) for Types 11, 111, and IV; in particular, one finds the following: For Type 11 there is a unique 
solution fii such that ^2(11) < hi < ^i(ii); there is a unique solution fm with fi(iii) < fm < f2(iii) for Type 111; and for 
Type IV there is a unique solution fiv such that ri(iv) < hv < ^2(iv)- 
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To summarize our findings here, we have found that the complete set T of stationary points of (20), including the 
singular stationary points, in the 6-dimensional, si, S2, s, a;i, a;2, a;-(phase) space Xdk, is the following 1-parameter 
family: 

J^ = J^{0 := = r2(jj,.s2(i) = rli^i),0,t^,^±rnj , (si(ii) = f'^n) , S2(ii) = ^2(„), 0, ^, ^, ^ ± jjn^ , 

(si(m) = ri(iij),S2(iii) = ?^2(iii)>0>'?>^.? ±mn) , (si(iv) = '^i(iv) > ■S2(iv) = f2(iy),0,^,^,^ ±r]ivj 
(si(i),S2(i),si = ri,^,^,^) , (si(ii),S2(ii),sii = rji,^,^,?) , 

(si(III)5S2(III)iSlII = ^^IIdCj^jC) > (si(IV))S2(IV))SlV = ?^rV>$)C)0 

(si(I),S2(I),Sl(I),^,^,0 ' (^l(I)i^2(I),S2(I),^,^,^) , (si(II),S2(II),Sl(II),$,^,^) , 

(si(II),S2(II),S2(II),^,^,0 ' (^1(III)> ^2(111), ' (^1(III)> ^2(111), S2(III),^,^,^) , 

(si(IV)>S2(IV))Sl(IV))^)^)0 ' (^1(IV))S2(IV)>S2(IV))^)C)0 : ^ € R} (55) 

5.2 Characterization of stationary points 

Owing to our analysis in Subsection 4.2, and the nature of the system of differential equations under investigation, it 
suffices to determine the type of the stationary points with respect to one of the invariant half-planes of the following 
forms, depending upon which stationary points are chosen for the pair of coaxial rings as specified by the four cases 
above: 

'^i (0 ■= {P = {si{i),S2(i),s,^,i,x) : s > 0, X e R} , 
'^ii (?) ■■={?= {si(u), 32(11), s,(.,^,x) : s > 0, a; e K} , 
Sin (0 := {P = (si(m),S2(iii),s,C,C,a;) : s > 0, a; e M} , 

^iv (0 := {P = (si(iv),S2(iv),s,?,?,a;) : s > 0, a; e M} (56) 

To do the Unear analysis, we have to compute the derivatives of the vector field in these half-planes at the stationary points 
Pe 5(0; namely, 

X'(P) := " ^ ' ^ M . (57) 
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We shall do this in rather complete detail for Type I (with a > 1), then in somewhat less detail for Type 11, and then just 
briefly sunmiarize the analysis for Types III and IV. Toward this end, we readily obtain the following formulas: 



(P) 



2{x-0 



9,$(P)=4rV / 



2 r 



I-/2 Kkfk 



7I-/2 Kkfk 



{x - - {r - h) (2r + fk) - 2rfk sin^ a 



cos 2ada 



{f — fk) + (a; — ^) + 4rffc sin^ (T 



5/2 



{r - fk) + 4rffc sin^ cr - 2 (a; - 



cos 2(7 da 



{r-fk) +{x-^) +4rffc singer 



5/2 



2 Kfef ^ I (2r^ — rfk —f^) ~ ~ 0^ + 2''?^fc sin^ cr cos 2(j — 6rfk sin^ crj rfcr 



2 ^7r/2 

a,*(P) = -6(a;-0E 

fc=i' 



(r — ffc)^ + (a; — + Arfk sin^ cr 
^^fc^^fe (?^fe — T cos 2(t) dcr 



5/2 



{r-fk) +4rffeSin^CT 



5/2 ■ 



(58) 



Now we fix a particular value of ^, which we denote as ^, and proceed to our Une of analysis of the stationary points in 
Si . First, we investigate the derivative at the fixed point pair P = pij_ := ^Si(i) = r^j^ , S2(i) = ^2(1) > 0> I i Vij ■ 
We compute from (58) that 

' 



(59) 



where 



1(1) 



1(1) + '/I 
^l(I) 



5/2 



+ 



2(1) 



2(1) 



5/2 



^3 



(^l(I)+^l)'^' (^2(I)+'?l) 



5/2 



(60) 



3x* {pi±) = TSttt?! 



f2 

^l(I) 



'2(1) 



(^1(1)+^?)'^' (^2(1)+'??) 



5/2 



where it should be noted that the partial derivatives with respect to s need to be calculated in a limiting sense as r J, 0. 
Whence, we immediately conclude that both pi^ and pi_ are saddle points, with stable and unstable manifolds and 
VF" in iSi (^), respectively, satisfying the following properties : 

W'{pi^) = (x|-\oo), 

and {pi+) at pi_^ is tangent to the line with nonpositive slope given as 

ii+ : S,* (pi+) s + [5^* (pi+) - 9,$ (pi+)] a; = 0, 
while ) at pi_ is tangent to the line of nonnegative slope of the form 

Li_ : 5,4- s + [a,* (pi_) - (pi_)] a; = 0. (63) 



(61) 



(62) 
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Moreover, it is easy to see from the form of the dynamical system (48) that VT" (pi^ ) \ {pi^ } ~ W" {pi_ ) \ {pi_ } :~ 
C\ := C\ {p\j_,p\_) is a smooth, convex curve joining pi_ to pi^, is symmetric with respect to the line x = S^, and is 
(asymptotically) tangent to Li^ and Lj at pij_ and pi , respectively (see Fig. 1). Thus, Ci together with the line segment 
{(0, x) : xj"-* < X < xj"''' } comprises a heteroclinic cycle that we denote as Zi, which is susceptible to perturbations that 
can generate chaotic dynamical regimes, such as described in ||4l, (|9], iflTl . Il25l . ||3T1 . BTl . \A2l . L53i . and [i6 1 1 . 

Next we consider the stationary points in !F of the form Q = (^si, §2, s, |, , where we have chosen a specific 
value ^ of ^ for our computations. As we are considering Case 1, si = Si(i) and §2 — 52(1)- Again we use formulas (57) 
and (58) to calculate the (linear) type of the stationary point via X'{Q). The following results may readily be verified: 



(Q) = 9.* (Q) ^ 0, 

2 p7r/2 



fc=i ■ 



—2' (64) 



1 .A /"^/2 (Sr^ - rffe(i) - f + 2rfj,(i) sin^ cr^ cos 2a - 6rfi,rT^ sin^ a 



da 



^fc(i)) + 4rffc(i) sin 



5/2 



Observe that, as was to be expected, these formulas are independent of ^. 

It is easy to see that dx^ (Q) is positive for all r > (and also all ^), and it can be shown by a routine but rather 
lengthy analysis that (Q) is positive at Q = gi := ^si(i), S2(i) , siii,i.,i^- Accordingly the eigenvalues of X'{qi) are 
±\/dx^ (qi) 9s (<7i), so qi is a saddle point. As for the singular stationary points Q = qi_ := (^si(i), S2(i), Si(i), ^, ^, 
and Q = qi^ := ^si(i) , S2(i) , S2(i) jiiiiCjy it is intuitively clear and follows from a straightforward analysis of the (48), 
in which we employ the results obtained in Subsection 4.2, that these points behave, except right at the points, just like 
centers of the Hamiltonian dynamical system. 

We may now, after some additional routine analysis of the system, piece together the results obtained on the stationary 
points and invariant curves in order to produce a complete qualitative picture of the nature of the phase portrait of (48) in 
any of the half-planes Si {£,), which is illustrated in Fig. 1. The main details of the phase planes includes the heteroclinic 
cycle Zi connecting the saddle points pij_ on the x-axis, which encloses the remaining three stationary points qi^, qi_, 
and qi and all of their associated stable and unstable manifolds. In particular, we see that 

[{qi}UW%qi)UWiqi)]n{r<h} 
forms a small homoclinic loop £i_ enclosing qi_ , and 

[{qi}UW%qi)UW"iqi)]n{r>h} 
forms a small homoclinic loop £i_^ enclosing gi^ . Moreover, it is not difficult to verify that 

[{qj}UW%qi)UW^{qi)]n{r>0} 

forms a figure eight curve, which we denote as (i, enclosing the singular centers qi and qi^. Owing to the symmetry of 
the system with respect to the line x — the homoclinic loops £i_ and and the curve (i are also symmetric under 
reflection in this line. In addition, the interiors of ii_ and are comprised of periodic orbits circling qi_ and qi^, 
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respectively, while points close to Ci and in its exterior belong to periodic trajectories that circle around it. Observe that 
this representation of the dynamics of (48) is quite consistent with what one might expect for a stationary, axisymmetric 
vortex breakdown flow of B-type. Thus Type I (when a > 1) we have, with our two coaxial vortex ring model, produced 
a two-parameter (a and k) family of axisymmetric B-type vortex breakdown flows. These features are illustrated in Fig. 
1. 

For Type 11, we find many similarities, but also some differences in the phase (half-) plane Su in comparison to 
the phase plane Si (5) for Type 1. Once again we have saddle points, pu_^ , on the x-axis determined from (54) with the 
obvious modifications, where their linear stability is determined via (60) mutatis mutandis. We also have the singular 
points, this time gii+ := (si(ii), S2(ii), Si(ii), I, and := (^si(ii), .§2(11), ,§2(11) , ■C, ^) . behaving like nonUnear 

centers. In addition, there is another stationary point gn := (si(ii),S2(ii),sii,|,^,^^, which lies between gn_ andgn+, 
and is a saddle point in virtue of (64) with the obvious modifications. There is also a heterocUnic cycle Zu (describing 
the bubble)cormecting pii_ and Pu_^, which is symmetric with respect to a; = ^, and encloses the stationary points qii_ 
and qu+ . The bubble shape produced for this type of fixed ring configuration tends to look less like the one expected, as 
compared to Type I for a > 1 sufficientiy large and Type IV as can be seen from Figs. 1, 2 and 4. 

Type III is almost completely analogous to Type II, where the only changes involve substituting si(iii) for si(ii) , S2(iii) 
for S2(ii), and sm for su. In particular, the phase portrait in iSm is quaUtatively the same as in Type II, consistent 
with the changes in the stationary points that we just indicated. However, we see from Fig. 3 that the bubble shape 
produced for this type of fixed coaxial vortex ring configuration tends to look less like the physical B-type flow than any 
of the other types, as can be confirmed by noting that the shape in Fig. 3 resembles a marble balanced on the axis of 
symmetry. Actually the shape for Type IV, as depicted in Fig. 4, seems to be most hke what one expects of a physical 
vortex breakdown streamline pattern of B-type. 

We note that in Figs. 2, 3 and 4, the small rectangles in the upper right-hand corners are microscopic depictions of the 
much smaller scale streamline topologies near the fixed coaxial vortex rings that are closest to the axis of symmetry. 
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Figure 1: Streamhnes for Type 1, a = 5, k= 1.5, and si ~ 0.06, S2 — 0.94, s ~ 0.24, x± ~ ±0.45. 
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We should point out that for different parameter values of a, k, the structure of the flow could be dramatically different 
as shown in Fig. Q for a Type I configuration for a substantially smaller than unity, where a hairpin vortex structure is 
present. We note that in this "singular" case, the dynamics is most certainly inequivalent to any of the instances illustrated 
in Figs. 1-4. More specifically, the behavior of the stable and unstable manifolds is very different from those in the other 
phase portraits: In this case there is a single homoclinic loop l\ beginning and ending at q\, which encloses the singular 
center . In contrast with the other behaviors illustrated above, however, the portion of W (gi) lying below q\\ actually 
goes off to X = cx), while the component of below q\\ tends Xa x — —oo, with both of these semi-infinite curves 

remaining above the cycle Z\, as depicted in Fig. 5. The analysis of this singular case - and the bifurcation behavior that 
transforms the dynamics in Fig. 5 to that in Fig. 1 - must wait for a later investigation, as we shall not attempt to treat this 
in any further detail in the sequel. 

We note, however, that hairpin vortex configurations have been observed computationally in the vortex flow in a 
cylinder with rotating covers Q. When the two covers are co-rotating, a part of the parameter plane (defined by the 
Reynolds number and the aspect ratio of the cylinder) may exist where a hairpin vortex is present in a steady flow. This 
region of the parameter plane grows to a significant size when the ratio of the angular velocities of the covers is increased. 



S 




X 



Figure 5: Streamlines for Type I, a = 0.1, k = 1.5, and si ~ 895, S2 — 6914, s ~ 3860, x± ~ ±24. 

It is eminently clear from our analysis in this section that Type I for a sufficiently large ambient swirling flow and 
Type IV , among all four types treated, produce the type of kinematics most consistent with the B-type breakdown be- 
havior encountered in experiments and obtained from numerical solutions of both the Euler and Navier-Stokes equations. 
Consequently, we shall only consider Type I with a > 5 in the sequel (keeping in mind the close analogs for Type IV) , 
and take advantage of this by reducing or dropping all of the unnecessary subscripts and superscripts in order to simplify 
our notation. We note only that most of the results that we shall obtain for this type of coaxial vortex ring positioning - 
including the existence of chaotic regimes for slightly oscillating, time-dependent kinematics - have direct parallels for 
the other types - and especially Type IV - as well. 
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5.3 The direct approach revisited 

As indicated above, from now on we shall consider Type I (with a large enough to preclude the hairpin vortex behavior 
shown in Fig. 5) almost exclusively. So in what follows, we employ simplified notation such as s\ and S2 for si(i) and 
S2(i) , respectively, s for s\, v for v\, and x± for x'^^ . Now we shall look once again at the kinematics from the perspective 
of the direct approach formulated in Subsection 3.2, but here for the case when the rings are located at the stationary points 
(si, S2, C) as described above. Since the rings are stationary, the potentially time-dependent Hamiltonian system (27), 
or (29)-(30), is actually autonomous. In particular, owing to our assumption that the x-component of the ambient velocity 
field is just the constant —a, the equations of motion are 

' ' COS lada 



(r — fk) + {x — ^) + 4rffc sin a 



2 3/2' 



^ f^^^ (f k — r cos2a) da 
x = -a + 2}2i^kfk — (65) 

From the dynamical equations (65) in comparison with the last two equations of (48), we see that the analysis of the 
stationary points and invariant manifolds for the passive formulation in this case (in Subsections 5.1 and 5.2) immediately 
carries over to the direct approach. More specifically, if the rings are initially at (si,^) and (s2,C)' respectively, then 
these points are singular centers, and the other stationary points are {0,x±), both of which are saddle points, and the 
saddle point (■5, 1^ In addition, the stable and unstable manifolds of the saddle points are just the same, with the obvious 
identifications, as those described above for p± and q, and give rise to the identical heterocUnic cycle - Z - and homoclinic 
loops £+ and £- forming the figure eight curve C,. Accordingly the phase portrait for (27) is indistinguishable - modulo 
the obvious change in labels - from that portrayed in Figs. 1 and 4 for Types I and IV. 

It is clear from our analysis so far that for fixed ring configurations, the dynamics of the passive fluid particles (i.e. 
the kinematics) is quite regular; in fact, the system is LA-integrable since it is Hamiltonian, autonomous, and has just one 
degree-of-freedom. Consequently, the Poincare map of first returns to the meridian half-plane Sj, which we introduced 
briefly at the end of Subsection 3.2, cannot admit chaotic regimes in the fixed ring case, and so must also be regular 
Therefore, it is manifest that if we are to generate chaotic kinematic regimes in our model - which are of particular 
interest in this investigation - the coaxial rings cannot both be stationary. In other words, a necessary condition for 
chaotic dynamics of the passive fluid particles is that the ring dynamics be nontrivial. We shall explore these matters in 
considerably greater depth in the next two sections. 



6 Kinematics for Slightly Varying Rings 

A number of studies taken together such as those of Breuer 111 51 . Br0ns et al. lfT6l[T7l[T8l . Gelfgat et al. Il23l . Hartnack 
et al. 123, Holmes |[291, Krause |I33] |34l EH, Lopez & Perry |[39J, Rusak et al. ED, Serre & Bontoux (|49l, S0rensen 
& Christiansen ||5T| , Sotiropoulos et al. ||52l[53 l, and Weimer ||60l, present compelling evidence, and in just a few cases 
proof, that the kinematics in B-type vortex breakdown exhibits chaotic regimes, which may include those of Shilnikov 
type. However, it appears that these chaotic effects are all traceable to small perturbations that break the axisymmetry of 
the configuration. Our model remains perfectly axisymmetric, so apparently the only way to generate chaotic kinematics 
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is to allow some non-stationarity for the system. Here we develop a framework for studying the kinematics of our model 
with small variations in the ring motion, and leave to the next section a more detailed analysis of the possible chaotic 
regimes generated by slightly varying rings. We begin by analyzing the dynamics of the two rings for small motions about 
the fixed points. 

6.1 Dynamics for varying rings 

To represent small motions of the rings about the fixed points, we set 

Sfe = Sfc + (^fc, x'fc = 1 + V'fe, (1 < A: < 2) (66) 

and note that it follows from the LA-integrability of (16), or equivalently (17)-(18), that the functions ^pk^ i^k in solutions 
of (16) are at least quasiperiodic, and are periodic if the initial conditions are chosen properly. We shall choose the 
perturbation functions so that |(^fe| and |V'fc| to be very small in a manner that shall be specified more precisely in what 
follows. The three-dimensional hyperplane 

(a) = {(si, S2, x\,X2) ■■ s\ + K,S2 = si + KS2 := a} (67) 

of the phase space of (16) is, as pointed out in (26), an invariant manifold of (16). To simplify matters, we shall only 
consider (varying) solutions of (16) that lie in (a), which implies that 

</92 = -«;~Vi- (68) 

This choice is also motivated by the fact that an examination of (46) readily reveals that the center surface (p) associ- 
ated to any stationary point p of the ring dynamics is contained in (a). 

If we choose the varying solution (66), with |v5fe(0)l , \ipkiO)\ (1 < < 2) sufficiently small, so that it is initially 
on W (p), then our varying rings orbit will be periodic. However, we do not know exactiy what (p) is at this point. 
Of course, returning once again to (46) and (47), it is straightforward to show that the linear approximation Wf^^ (p) of 
W (p) is defined by the homogeneous Unear equations 

Sl -I- KS2 = 0, 

Axi - Bx2 = 0, (69) 

where 

A = A (si, S2, 1; a, := [v^ 4- 5^, $2 (p) (p) - nd^.^i (p)]} , 

B = B (si,S2, := $2 (p) (p) - Kds.^i (p)] , (70) 

and 

y = V5xi*2 (p) (p) - (p) + K (p) - ds,^2 (P))]. (71) 

Observe here that we can use the above to approximate the frequency uj of periodic solutions near the center in (p); 
namely 

w ~ z/. (72) 
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Moreover, if we select the initial point to be on W{i^ (p) and very near to, but distinct from the stationary point p, we have 
the following additional useful (but only approximate) relationship between the perturbing functions - defined in (67) - 
for the coordinates xi and X2'. 

^2 ^ ^^1. (73) 

All of the above suggests a simple method for locating a point (distinct from p) of {p): Choose small values for 
(fii (0) and (p2 (0) that satisfy (68), and select ipi (0) to be small and nonzero. Then, if we select a very small values of all 
of ipi (0) , (f2 (0) and -01 (0), we are certain to find the value of ■02 (0) satisfying (si(0), S2(0), xi (0), 2:2(0)) G {p) in 
a very small interval centered at (A/B) (0); for example, by using a simple bisection method. Naturally, once we have 
found such a solution, we are guaranteed that it is periodic and that its period is approximately equal to ^-k jv. 

For a more precise specification of the size of the perturbations around the fixed points, we turn to the following result, 
which is illustrated in Fig. 6. 

Lemma 1. The varying solutions of (16) sufficiently close to a point on the center manifold W (p), which is in a 
sufficiently small neighborhood of the fixed point p — ^§1, §2, i, can be expressed in terms the following expansions 
in terms of a small parameter fi: 

si — si [t] fi) ~ si + fismi't + O (/i^) , S2 = S2 — (m/k) sinj/t + O (/i^) , 

xi = xi {t; ^) = ^ + /i cos i^t + (fi^) , X2^i + f-L {A/B) cos i^t + O {i?) , (74) 

where the order relations are uniform in t when < e* :~ (V^) min{s — Si, S2 — s}, with si, S2,and s as defined in 
our discussion of (48) and (65) in the preceding section. This restriction guarantees - for one thing - that the orbits of the 
rings do not intersect one another 



Proof: This follows directly from the LA-integrability of the system, the properties of W{[^ (p) in relation to W {p), 
and the description of the stationary points and phase space structure for the corresponding dynamical systems delineated 
above. ■ 
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Figure 6: Slightly varying ring dynamics and stagnation point inside the bubble for a = 5, ^ — 0.01 (left) and a — 
20, /i = 4 X 10^'^ (right). The unperturbed heteroclinic orbit is also shown. 
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We note that in Fig. 6 the small rectangular boxes in the upper left hand comers depict the much smaller scale elliptic 
orbits of the slightly oscillating rings in neighborhood of the upper fixed points, while the small rectangular boxes to the 
right within the figures provide a microscopic view of the dynamics near the lower fixed points, which are extremely close 
to the axis of symmetry. 



6.2 Kinematics for oscillating rings 

Returning to (29)-(30), it now follows directly from Lemma 1 and a straightforward calculation, left to the reader, that we 
can expand (30) in powers of fj, as described in our next lemma. 

Lemma 2. For the varying ring dynamics as described in Lemma 1, we can expand the Hamiltonian function of (29) - 
representing the direct kinematics formulation - as 



H (s, X, t) = Hq (s, x) + iiH\ {s, x,t) + (/i^) , 



uniformly in t when < e*, where 



''/^ cos2ada 



/ ^ . - r' cos: 

Ho (s, x) := as - 4r > Kkrk / — r- 



1/2 
k 



and 



Tv/2 



ni{s,x,t) :=2rsmvty2{-'^f [ 

fc=i -^0 

4r (^x — cosvt J 



r{r — fk) + {x — ^ + 2rfk sin^ f 



cos 2crda 



7r/2 



- A 3/2 

fi KAr2 \ „ , 
H T-^TTK ? cos 2<jda. 



- + 



(75) 



(76) 



(77) 



With Lennmas 1 and 2 now at our disposal, we have the elements necessary to delve more deeply into the behavior of the 
kinematics - as characterized by the time-dependent Hamiltonian system (29)-(30). 

6.2.1 Oscillation of stagnation points 

When the rings are stationary, we showed in Subsection 5.1 that the stagnation points of the advected flow on the a;-axis 
are fixed at the points x± determined by (51) and (52). In contrast, when the rings are oscillating in accord with Lennma 
1, these stagnation points will also vary slightly in a quasiperiodic (or periodic if the rings are properly tuned) manner. We 
denote the time- varying positions of the stagnation points by x± {t; fi), noting that these are independent of t when fi = Q 
and given as 

x±{t;Q) = x±. (78) 
These points are determined in a manner analogous to x±; namely, they are the solutions in x to 



^ {x, t; /i; k) = —a + n 



Krl (t; ii) 



Si {t; ji) + {x - xi{t; ii)Y 



1 3/2 



S2 {t; jj) + {x - X2{t; ij)y 



3/2 



0, (79) 
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where Sk {t; /i) and Xk{t; fi), 1 < k < 2, are as defined in Lemma 1 . It is straightforward to show, just as for the stationary 
case described by (51) and (52), that, under the restrictions on fi specified in Lemma 1, (79) has precisely two solutions 
satisfying 

x_{t;fi) <i<x+{t;fi) (80) 

for all < G M, where {t; ^) and x^ {t; /i) are only approximately symmetric with respect to ^, rather than exactly 
so as in the stationary ring case. The oscillation of the stagnation points is illustrated in Fig. 7 for several parameter 
values. Naturally, we would expect this symmetry breaking to have consequences regarding the nature of the time- varying 
vortex breakdown bubble. As we shall see, this loss of axisymmetry can indeed have very dramatic consequences on the 
kinematics. 
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Figure 7: Oscillation of the stagnation point along the axis of symmetry (s = 0) for a — 5, fj, — 0.01 (left) and 

a = 20, = 4 X 10-^ (right). 



6.3 Natural Poincare maps for kinematics 

Let us return to the meridional first return Poincare map introduced at the end of Subsection 3.2, which provides a very 
natural and useful means for (discretely) characterizing the dynamics of passive fluid particles induced by the small 
oscillations in the positions of the coaxial vortex rings described in (74). For concreteness and simplicity, we shaU assume 
in the sequel that the rotating component of the superimposed swirling flow has the (constant) form 

= n, (81) 

with rt a fixed positive parameter; in other words, 6i = 62 = in (14). In this case, the Poincare map has a particularly 
simple definition; namely, 

77 = 77 (•, •; a, k, fi) = Uf, : Sj Sj, 

where 

(so, xo) (s {2n/n) , x {2n/D.)) , (82) 

and (s (t) , X (t)) is the unique solution of (27) satisfying (s (0) , x (0)) = (sq, xq). Note that although the Poincare 
map depends on all of the parameters a, fi, k, and fi (and incidentally on x), we have emphasized its dependence on the 



27 



parameter ji connected with the oscillation of the vortex rings by singling it out via a subscript. With this notation, //q 
corresponds to the fixed ring case associated to TLq in (75). 

As we shall see in our analysis of chaotic motions in the next section using Melnikov's method, it useful to list a 
number of elementary properties of the Poincare map iJ, all of which can be readily deduced from its definition in terms 
of the flow of the Hamiltonian system (29)-(30), which is equivalent to (27): 

(PMl) The map : ^ is a (real) analytic diffeomorphism (actually a symplectomorphism), which is analytic in 
both its variables and its parameters when they are restricted to the range of values delineated above. 

(PM2) iJo has hyperbolic fixed points at (0, x+) and (0, X-) connected by a heteroclinic cycle Z comprised of the line 
segment joining these points and the curve C as described in Subsection 5.2. Moreover, it has the following fixed 
points: (si, and (s2, C), which are centers of the symplectomorphism 77o, along with the saddle point (s, ^) . 

(PM3) If the dynamics of the rings is (nontrivially) periodic of period 2-11/10 and cj/il = n G N, the set of natural numbers, 
both (0, x+) and (0, xJ) are hyperbolic fixed points of whenever < e*. On the other hand, if fi/oj = m, a 
positive integer greater than one, both (0, and (0, x^) are hyperbolic periodic points of I/^i of (least) period m. 

(PM4) In any case, if /i satisfies < e for a sufficiently small < e < e* it follows from the axisymmetry and the 
standard results on the persistence of hyperbolic fixed points (see e.g. ll25l 1611 ) that 11^ has a pair of hyperbolic 
fixed points p+ (/i) := (0,i+(/i)) andp_ (/i) (0,2;_(/x)) such that: (i) £±(0) = x±; (ii) ai-(/i) < for all 

and (iii) the stable and unstable manifolds of these fixed points satisfy M^'' (p+ (/i)) ~ {(0, x) : x^{fi) < x}, and 
(/i)) = {(0, x) : X < Moreover, 77^ has a hyperbolic fixed point q (/i) := (s(/i), satisfying 



Additional insight into the nature of the map 77^ can be obtained from the results of several numerical experiments 
presented in Figs. 8 and 9 in Section 8. 

6.3.1 The heteroclinic cycle 

As we shall see, the heteroclinic cycle associated to the fixed ring kinematics plays a key role in the identification of 
chaotic kinematic regimes for slightly oscillating ring dynamics. Accordingly a more detailed analysis is in order. To 
this end, we recall our definition of the heteroclinic cycle, which joins the leading and trailing points of the bubble, and 
encloses all other stationary points along with their stable and unstable manifolds defining the homoclinic orbits £± and 
the smaller heteroclinic cycle (; namely 



where C± is the closed axial line segment in f) defined (for a fixed ^ = ^ as C± := {(0, x) : X- < x < x+}, and C is 
the convex, ^-symmetric curve beginning at (0, x+) and ending at (0, X-) as described in Subsection 5.2. 

In virtue of (65) and (PM2), we know that Z is a heteroclinic cycle for both the continuous , LA-integrable, Hamilto- 
nian dynamical system 




Z := C±\JC, 



(83) 



s = dxUois.x) , X = -dsHo{s,x) , 



(84) 
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where the Hamiltonian function is defined as in (76), and the discrete, symplectic (integrable) dynamical system defined 
by the map 

no-.Sj^S), (85) 

in accordance with (82). To find the heterochnic orbit in £±, it follows from (65) that we must solve the differential 
equation 



X = —a + IT < 



3/2 



ri + 



3/2 



(86) 



which is no easy matter. Fortunately, we do not need a complete closed form solution. To illustrate this, let the heterochnic 
orbit satisfying a; (0) = ^ be denoted a.s^i{t) : = {si{t),xi{t)). Then, it follows directly from its definition that si{t) = 
for all t e M, and 

xi{t) =x-+0 (e^*) and xi{t) = O (e^*) as t -oo, (87) 



and 



xi{t) =x+ + {e-"*) and xi{t) = O (e"^*) as t ^ oo, 



(88) 



where x is a positive constant. 

On the other hand, it is clear from (76) that C is the solution curve of the equation 

2 



ar 



2 ^7r/2 

fe=i -^0 



cos 2ada 
A 1/2 



= 0, 



(89) 



which is obtained by factoring r out of the energy curve Hq {s, x) = 0. It is possible through a rather deeper analysis of 
(89) to extract many details concerning the form of C, but it turns out that we just need to use some simple asymptotic 
representations for the heterochnic orbit := {su{t),Xu{t)) lying in C and satisfying the initial condition s„(0) > 

0, Xu{0) = ^, along with a few elementary symmetry properties of this orbit. These could be determined by solving the 
differential equation obtained from (67) and (89) - a rather formidable task. Fortunately, the asymptotic and symmetry 
properties of that we require are rather easy to deduce. Owing to the definition of the heteroclinic orbit, it is easy to 
verify that 

Xu{t) =x+ + (e^**) and i;„(t) = O [e"*) as t ->■ -oo, (90) 



and 



with 



and 



Xu{t) =X-+0 (e""^*) and = O (e"''*) as t ^ oo, 
s^{t) = O (e^*) and s„(i) = O (e''*) as t -oo, 
s„(t) = O (e-^*) and s„(t) = O (e"^*) as t ^ oo. 



(91) 



(92) 



(93) 



where the positive constant x may have to be adjusted in order to satisfy all six expressions (88)-(89) and (90)-(93). As 
for the symmetry properties, it is clear from the relevant definitions that the following obtain: 



Sui-t) = Suit) > and Xui-t) - ^ = -(a;„(t) - 0, (with - < for t > 0) 



(94) 
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and 

Su{-t) = -Su{t), (with s„(i) < for i > 0) and x„(-t) = x„(t) < (95) 

for all — cxo < t < oo. 

7 Chaotic Kinematics 

We shall demonstrate in this section that small oscillations in the coaxial vortex ring motions lead to chaotic dynamical 
regimes for the advected fluid particles: in short, small oscillations generate chaotic kinematics. This shall be accom- 
plished by using Melnikov's method to demonstrate that slight motions of the rings can lead to transverse intersections 
in the heteroclinic cycle Z, which generate chaotic regimes in virtue of well known results such as presented in Wiggins 
||6T1 . We note that the axisymmetry of the dynamics and kinematics precludes such intersections in the C± portion of the 
cycle, so the only possible such intersections must occur between the unstable manifold W'^{p+ (/i)) and stable manifold 
W{p- (/J.)), which comprise C in the stationary ring configuration for /.i = 0. 

The chaotic behavior can be summarized in the following result, which also can be proved by making a fairly straight- 
forward - but far from obvious - modification of the three coaxial ring analysis in Bagrets & Bagrets lIU . 

Theorem 1. For each < < e, the map 

defined by (82) has a heteroclinic cycle (with Zq ^ Z) having transverse intersections, which implies the existence of 
chaotic orbits. 

Proof: The key to the proof is an analysis of the zeros of the Melnikov function (cf H, ID, |[TT|, ED, ED, ED, ED, 
andlMl) 

/OO 
{no,ni}^{(t>u{t),t + T)dt, (96) 
-oo 

where the standard Poisson bracket is defined as usual as 

{7^o,Wi}o d^UadsHi - dsHod^Hi. (97) 
It is clear from (29) that the Melnikov function may be written as 

/oo 
[su{t)dsHi + r) + xu{t)d..ni + ^)] dt, (98) 

-oo 

and this can readily be shown to be well defined owing to the convergent nature of the integral, which follows directly 
from (84) and (90)-(93). 

Now we compute from (77) that 

9s '^^1 (0n it),^ + t) — 01 (t) (sin i^t cos vt + cos vt sin vt) + 62 (t) (cos i^t cos — sin vt sin i^t) , (99) 

and 

dxTi-i {4>u{t),t + t) — 2ru {t) [63 (t) (sin j/f cos i^t + cosj^tsini/r) + 64 (i) (cos z^i cos i^t — sinz/i sinz/r)] , (100) 
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where 



ei it) (-1)' ^k' C'^ K"^ (Mt)) \r-\t) (x4t) - i 



62 (t) 



{ru{t)-fk) +[xu{t)-n + 



ru{t)fk {2 -cos2a)]-2fkCOs2a (r„(i) - f^)^ - 2r„(i)rfe sin^ cr | 



cos 2ijda, 



(101) 



7r/2 



(i ) - ?y - (r„ (t ) - f 1 ) (3r„ (t ) + fi ) ■ 



(ruit) - h) (3r„(t) + 7=2) - 2ru{t)f2 sin^ ct] } cos2CTdCT, 



(102) 



and 



fc=i 



[2r„(t)ffe sin^ cr- 



64 {t) := 2 



7r/2 



{ru{t) ~ fk) {ru{t) + 2fk) - [Xu{t) - i 



cos 2ada, 



(103) 



{ru{t) -h f ~2 (^Xu{t) -t)^ + 4r„(t)fi sin^ a 



+ 



iru{t) -Hf -2 (xu{t) - ' + 4r„(t)r2 sin^ a 



cos 2ad<7. 



(104) 



Observe that it follows directly from (94) and (95) that 6i and 83 are even functions of t while 62 and 84 are odd 
functions. 

Our analysis so far, together with the identification of odd and even functions, leads immediately to the following 
simplification of the Melnikov function: 



ajl (r) = 2 cos I'T / {[s„(t)ei(t) + 2r„(i)i„(t)e3(i)]sinzyi- 
^0 

[Su{t)e2{t) + 2r„(t)i„(t)e4 (t)] cosiyt} dt. 
Whereupon, a sti-aightforward but rather tedious analysis of the integral shows that 

971 (t) = Ccosz/r, 



(105) 



(106) 



where C > 0. Consequently, 9Jl (t) has simple zeros, which proves that there is a transverse heteroclinic orbit in the cycle 
Z^, and this implies the existence of chaotic dynamics for the discrete system generated by 77^ for /i ^ 0. Thus the proof 
is complete. ■ 

A few remarks are in order concerning the nature of the chaotic kinematics described in Theorem 1 . We have shown 
that small oscillations of the rings leads to a transverse intersection in the large heteroclinic cycle Z, but one can also 
prove that such motions induce transverse intersections in the homoclinic loops l± comprising the figure eight curve 
C(cf. lIU). Consequently, it is not surprising that there is an accumulation of chaotic streamlines around the stationary 
saddle point (s, as seen in several of our figures. Of course, the chaotic behavior we have proved cannot - in virtue of 
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axisymmetry - include a regime of the type discovered by Shilnikov ll50l and which has now become a standard fixture 
in modern dynamical systems theory (c/. 1251 [311 1611 ). However, a close look at the chaos resulting from a transverse 
intersection of the stable and unstable manifolds comprising the upper branch C of the heteroclinic cycle Z, reveals certain 
Shilnikov-like features of the dynamics of the passive fluid particles. In particular, denumerably many orbits through 
points arbitrarily close to the axis of symmetry C± (which contains the unstable manifold of the trailing point 
described in (PM1)-(PM4)) of the bubble, accumulate chaotically around the trailing and leading points of the bubble, 
which is a characteristic and fundamental property of the dynamics of the three-dimensional Shilnikov model. We note 
that dynamical systems considerations and several experimental and numerical studies of fully three-dimensional B-type 
vortex breakdown regimes subject to non-axisymmetric perturbations indicate that Shilnikov chaos is quite common in 
the streamline patterns within the bubbles lITSl l29l l52l [53l . 



8 Illustrative Examples 

In this section we study several cases that demonstrate the kinematic behavior of our model (for the Type I positioning 
of the rings) via the Poincare map 77^, and provide insights into the ways in which it depends on the various parameters, 
especially with regard to the onset of chaotic regimes. We shall for each of the cases take x = Ij 000, which is tantamount 
to prescribing the core radius of the rings to be 0(0.001). To give a fairly representative range of possibihties, we shall in 
particular compute and present pictorial representations of the dynamics of 77^ for the following cases: 

Case 1: a = 5, 1) = 1/ ~ 7.96, k = 1.5; (a): a* = 0; (b): fi = 0.001; (c): fi = 0.01. 

Case 2: a = 20, = ~ 199.7, k = 1.5: (a): fi ^ 0\ (b): ^ = 4 x IQ-^; (c): fi = A x IQ-^ 
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Figure 9: Poincare section for a = 20, k — 1.5, ft 



199.7 for a) /X = 0, b) ^ = 4 X lO"'"^ and c) = 4 x 10" 



The Hamilton equations are solved using a 4th order Runge-Kutta scheme with a time step of the order of 10^^. 
From the result of the previous section, one expect the homoclinic and heterocUnic orbit to break after a small perturbation 
is added to the system as seen in Figures [8] and |9] We see from these figures that as the size of the perturbation increases, 
transverse homocUnic intersections interior to the bubble and transverse heteroclinic cycle behavior of the bubble bound- 
ary becomes more pronounced, so more regular tori are destroyed, leading to larger chaotic regions inside the bubble - 
indicated by the characteristic splattering, along with islands in a sea of chaos. Note in particular how the heteroclinic 
orbit delineating the upper boundary of the bubble appears to create a chaotic layer surrounding the bubble, and that the 
thickness of this layer increases with the size of the amplitude of the ring oscillations. One also should notice the escape 
of passive particles at the tails of the bubbles indicated by the arrows, which presumably occurs in devil's staircase fashion 
(cf. 11811281153]) 
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9 Epilogue 



We have provided a rigorous demonstration that a pair of coaxial vortex rings immersed in a swirhng ideal fluid flow appear 
to generate all of the kinematic behavior associated with vortex breakdown of B-type - including chaotic streamlines with 
some Shilnikov-like features when the rings undergo small oscillations. Naturally, although our results may turn out to 
be very useful in studying vortex breakdown, we freely admit that we have by no means proved that multiple coaxial 
vortex rings are the fundamental underlying mechanism in vortex breakdown phenomena. But even if the multiple rings 
do turn out to be the ultimate engine of vortex breakdown, one may well ask - as we have asked ourselves - how are these 
rings formed in swirling flows? The answer to this question might turn out to be quite elusive; nevertheless, we shall next 
endeavor to posit a plausible explanation that combines some mathematics and basic fluid mechanics. 

To see how the two slender coaxial vortex ring structure (described in Types I - IV, and perhaps most aptly for Type IV 
configurations) might develop, we start with only the swirling flow having its vorticity (vector) along the axis of symmetry, 
which we have taken to be the x-axis in our analysis. At the instant when a breakdown bubble is formed, the projected 
flow on any meridian half-plane must have at least one stationary point within the bubble owing to the Poincare-Hopf 
index theorem ||3l |25l [3T1 |6T| . As the flow structure of the bubble persists even in non-Hamiltonian real (viscous and 
compressible) fluid dynamics, this observability of the B-type breakdown phenomenon is suggestive of structural stability 
for the bubble formation process, which requires that the initial stationary point must be a saddle, associated, in virtue of 
the axisymmetry, with the intersection of any meridian half -plane with a circular streamline about the x-axis corresponding 
to a hyperbolic periodic orbit of the flow. But then, once again by the Poincare-Hopf index theorem, at least two centers 
are required to balance the saddle point. Inasmuch as the process of bubble formation is essentially instantaneous, one 
would expect the necessary two centers to be singular, which is consistent with the appearance of a pair of coaxial vortex 
rings. There is also a vorticity conservation element in this process in accord with Kelvin's theorem lfT9l[36ll44ll47ll57]| . 
which may be explained along the following fines: Tfie appearance of tfie bubble alters the velocity field in a way that 
introduces a vorticity field transverse to the initial field along the x-axis. The induced transverse vorticity component 
that results from the bubble formation must be balanced by a structure within the bubble such as a coaxial vortex ring of 
nonzero strength, and as we have seen there must be two rings in order to satisfy the index theorem. 

Some of the flow structures found here are richer than those found computationally for the flow in a cylinder with 
rotating covers ll5l [T6l . The general appearance of a hyperbolic closed streamline found in the present study gives rise to 
structures which, borrowing a notion from lfT6l . are "bubbles with an inner structure" as shown in Figs. 1-4. These are 
only rarely found in the flows in confined cylinders. However, as the hyperbolic closed streamline occurs quite close to 
the axis, the overall flow patterns do not deviate much from those observed in real flows, even if from a strict topological 
perspective they are quite different. 

Although our scenario for two coaxial vortex rings driving the B-type vortex breakdown is reasonably consistent with 
the governing physical laws and observed behaviors, we have no illusions about it being indisputable. Be that as it may, 
we have shown that such a model does provide a rather complete description of B-type breakdown flows. Accordingly our 
approach has the potential for providing an excellent paradigm and vehicle for developing techniques and instrumentalities 
for ameliorating, focusing and controlling B-type vortex breakdown phenomena, and we hope to investigate some possible 
applications in these directions in the not too distant future. 

During the course of the analysis of our two ring model, we were struck by the strong similarities between the kine- 
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matics induced in a meridian half-plane and the phase plane structure of the kinematics generated by a pair of point 
vortices in a half-plane. This is hardly surprising in view of the strong connections between coaxial ring dynamics and 
point vortex dynamics estabhshed in such investigations as ifTOll . |fT3]| . and lfT4l . Moreover, it has been proved ||9l that the 
motion of a pair of point vortices in a half-plane is capable of generating chaotic behavior in the advection of passive fluid 
particles. Given the simplicity of the kinematics generated by a pair of point vortices in a half-plane compared to that 
of a pair of coaxial vortex rings, which is essentially on the order of logarithmic versus elliptic integral Hamiltonians, it 
seems natural to investigate the possibility of the former - perhaps immersed in a parallel flow - producing the complete 
spectrum of B-type vortex breakdown kinematics. The practicability of employing the two point vortex model as a more 
tractable paradigm for vortex breakdown phenomena of B-type is quite an attractive option, which definitely calls for 
further exploration. 

Speaking of natural problems and questions for further study, three that we have avoided in this paper come to mind. 
First, ambient swirling flows with translational and rotational components that are dependent on the radial distance from 
the axis of symmetry are certainly more representative of the behavior of real fluids, where one expects to see velocity and 
rotation profiles that tend to decrease with distance from the axis. Models incorporating such radial dependence should be 
investigated, and one expects to be able to obtain even more realistic B-type breakdown behavior in such cases. Secondly, 
although we pointed out the possibility of rather unexpected hairpin vortex + bubble structures, seeming to appear from 
"out of the blue", associated with Type I for small translational velocities in the ambient swirling flow, we made no 
attempt in this paper to analyze them. There is clearly some very interesting - and challenging - bifurcation behavior 
associated with these unusual structures warranting further investigation. In particular, how does the hairpin + bubble 
geometry bifurcate into a single bubble as the ambient translational velocity increases? Lastly, it might be interesting to 
study analogous vortex breakdown models for bubble structures based upon three or more slender coaxial vortex rings 
immersed in swirling flows, again in aid of obtaining even more realistic B-type vortex breakdown structures. 

Perfect axisymmetry is, of course, as much an abstract ideal as is an ideal fluid, and as we have already remarked, there 
have been a variety of non-axisymmetric behaviors observed in experiments and numerical studies, which are precluded 
by the enforced axisymmetry of our model. It is therefore quite natural to enquire whether or not a full three-dimensional 
version of our coaxial vortex ring model that admits small symmetry breaking perturbations is capable of producing all of 
these non-axisymmetric kinematic effects? We plan to undertake such an investigation in the near future, fully expecting 
to see all of the non-axisymmetric dynamics, including the emergence of vortex breakdown of S-type at the traihng edge 
of the bubble - which may signal transition to turbulent flow regimes, and Shilnikov chaos. 

Finally, there are those who might simply say that the ideal fluid flow context (enabling us to reap the considerable 
benefits of a Hamiltonian structure for the kinematics) of our B-type breakdown model disqualifies it from serious consid- 
eration of such phenomena in real fluids. Certainly there is considerable truth to this with regard to some flow quantities; 
however, extensive experimental and numerical studies of B-type vortex breakdown in both real, high Reynolds number 
flows and ideal flows have shown that both B-type and S-type vortex breakdown behavior tends to be qualitatively the 
same in all of these cases, and even quite close quantitatively speaking over time intervals of several minutes. Thus the 
ideal tells us a great deal - especially dynamically - about the real when it comes to vortex breakdown. 
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